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Abstract. - In the dynamics of the undamped Frenkel-Kontorova model with kinetic terms, 
we find a transition between two regimes, a floating incommensurate and a pinned incommen- 
surate phase. This behavior is compared to the static version of the model. While in the static 
case the two regimes are separated by a single transition (the Aubry transition), in the dynam- 
ical case the transition is characterized by a critical region, in which different phenomena take 
place at different times. In this paper, the generalized angular momentum we have previously 
introduced, and the dynamical modulation function are used to begin a characterization of this 
critical region. We further elucidate the relation between these two quantities, and present 
preliminary results about the order of the dynamical transition. 

s 

Introduction. The Frenkel-Kontorova (FK) model Jlj] describes the interaction of a 
harmonic chain of atoms with a rigid substrate with period incommensurate to the lattice 
parameter of the chain. Its generality makes it a powerful model to investigate many different 
physical systems (^,[|, and in particular microscopic friction between contacting surfaces ^-^|. 
The static version of the model is characterized by the Aubry transition from a floating 
to a pinned state, for a critical value A c of the substrate modulation potential. Using the 
undamped dynamical version of this model, we addressed in previous papers the topic of 
"dissipation" (in the sense of transfer of energy from the center of mass to phonon modes) 
in incommensurate structures: we have studied the mechanism (parametric resonances) that 
governs the onset of sliding friction ||, and the conditions under which a new conserved 
quantity can be denned, which can be seen as a Generalized Angular Momentum (GAM) in 
the complex plane [||. 

In this work, we present new results, showing that, in the dynamics, a floating-to-pinned 
transition, analogous to the static Aubry transition, is found for all values of the potential 
A < A c . The transition is characterized by a region of critical times, with a remarkably 
complex behavior. After describing the FK model, we show that the dynamical modulation 
function undergoes at a critical time t c ± the same breaking of analiticity that, in the static 
case, occurs at A c . A second critical time t C 2 > t c i, at which the GAM conservation stops, 
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Fig. 1 - Left: the modulation function in the static FK model, (a): for A < A c the function is analytic; 

(b) for A > Ac it becomes non-analytic. Right: two snapshots of the dynamical modulation function. 

(c) TV = 233, A = 0.05, t = 300 < t& ; (d) N = 233, A = 0.05, t = 820 > t cl . Even if A < A c , the 
modulation function becomes non-analytic. 



was identified in [pj; here, we show explicitly the connection between these two quantities, by 
showing analytically that the analiticity of the modulation function implies conservation of 
the GAM. We also present some initial results about the order of this dynamical transition. 

The static model. The FK model describes a chain of N atoms in contact with a 
substrate, described by a rigid periodic potential. The intrachain interaction is modelled by 
a first neighbor harmonic potential, and the total potential energy reads: 



V = E 



r [u n +i ~ u n - I) +— sin ) 

2 Z7T Til 



(1) 



where u n are the lattice positions, I is the average atomic distance, and A is the strength of 
the substrate potential rescaled to the spring constant. The ratio between / and the mod- 
ulation potential period m is chosen to be an irrational number (in our case, m = 1 and 
I = t = (v5 + l)/2, the golden mean), in order to simulate an incommensurate (IC) modu- 
lated structure. Numerically, we use a finite chain of length Nl and use ratio of consecutive 
Fibonacci numbers to approximate incommensurability Q . The Aubry transition can be mon- 
itored by the modulation function (see e.g. ||). It is defined as: f(nl mod m) — u n — nl — Q, 
where Q = -h ^ n u n is the center of mass (CM) of the chain, and contains information about 
the groundstate positions of the chain particles with respect to the substrate. When A > A c , 
it becomes non-analytic, as shown in the left side of fig. |l|. Another feature of the FK model 
is the existence of a zero frequency phason mode below A c that moves to finite frequency at 
the Aubry transition. If the interaction with the substrate is weak (A << A c ), deviations 
from equidistant spacing / in the groundstate are modulated with the substrate modulation 
wave- vector q = 2irl, due to the frozen- in phonon uj q . For values of the coupling A < A c , 
the Fourier transform of the displacement related to q has thus the strongest amplitude, and 
amplitudes of higher harmonics nq scale with A'™'. It is therefore natural to use |n| instead of 
k and relabel the modes (see || for the details of the relabeling procedure). In the left side of 
fig. |2|, it can be seen that, both with the usual labels (panel (a)) and with the relabeling (panel 
(b)), the effect of the modulation is lost when A > A c , which for this model is A c = 0.154.... 
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The dynamical model. - In the dynamical version of the model, a kinetic term is added 
to the potential energy; in our case, it corresponds to giving an initial velocity to the particles 
of the chain. The Hamiltonian now reads: 

h = v + jc = v + J2y ( 2 ) 

n 

where the p n are the particle momenta. The CM velocity is defined as: P = -k ^2 n Pn- The 
dynamics of the system can be studied by following the time evolution of the CM of the 
chain Q, and of the particles deviations x n = u n — nl — Q. As initial conditions we choose 
Pn(t = 0) = Pq and x n (t = 0) corresponding to the groundstate. The equations of motion for 
the deviations are numerically solved with a Runge-Kutta-Fehlberg algorithm. Convergence 
of the results is found for a tolerance of 10 -9 , ensuring total energy conservation up to 10 -8 . 
As we saw in the previous section, phonon amplitudes play a crucial role in the statics. 
In the dynamics, their behavior can be used to elucidate the time evolution of the system. 
The inherent nonlinear coupling of the CM motion to the phonons of the chain leads to an 
irreversible decay of the CM velocity, that in an Hamiltonian picture can be identified with 
the onset of sliding friction. Transfer of kinetic energy from the CM to the internal vibrations 
occurs via a complex sequence of parametric resonances mediated by ordinary resonances with 
phonons related to the modulation potential 0] . These excitations can further combine and 
give rise to the appearance of Umklapp terms. This phenomenon can be monitored by the 
breakdown of the conservation of the GAM, as derived in ||: 

P<t> = ^ n 9 a -nga„ g + — Q = L + — Q (3) 

n 

where a± nq are the phonon amplitudes, obtained by Fourier transforming at each timestep 
the expression for the chain distorsions x n . 

Modulation function and phonon amplitudes. - We have shown in the previous section 
how the behavior of the modulation function changes when going through the transition at 
A c . If we now turn to the dynamical FK model with initial velocity, it is very intriguing to 
see that, for all values A < A c , the same behavior is found, after a critical time t c \. 

This can be shown by extending the concept of modulation function to the dynamics. 
Since the value q = 2wl is fixed, we do it in a straightforward way: ft(nl mod 1) = u n (t) — 
nl — Q(t) = x n (t). The results are shown in panels (c) and (d) of fig. |l|. The similarity 
with the static case is apparent. The time t c \ can be accurately determined from the largest 
discontinuity of the modulation function (see fig. |], lower panels) . Another signature of the 
dynamical transition is that it is accompanied by a breakdown of the conservation of the GAM. 
This is shown in the middle panels of fig. [}| It is however remarkable that this breakdown 
always occurs at a later time i C 2 > t c \. Fig. |J shows that this remains true for all system 
sizes considered. Qualitatively, one can understand this because the non-analiticity of the 
modulation function is caused by the excitation of so many phonon modes as to render the 
Fourier series representing it not absolutely convergent. The conservation of the GAM only 
breaks down when this excitation reaches the zone boundary, generating Umklapp terms ||. 
This is shown more formally in the next section. This dynamical transition can also be seen 
in the qualitative difference in the Fourier amplitudes shown in panels (c) and (d) of fig. |2|. 

The GAM and the modulation function. - We show that analiticity of the modulation 
function implies conservation of the GAM. The modulation function can be written as: 

f{z) = Y^a mq e a ™ z (4) 
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Fig. 2 - (a) Phonon amplitudes in the static FK model, for the cases N = 233, A = 0.05 < A c , and 
A = 0.16 > A c . The modes are labeled in the usual way. The first two modulation wave- vectors 
k — q,2q are shown, (b) Same, but the relabeling k — > n is used to make the scaling behavior 
apparent, (c) Phonon amplitudes in the dynamical case, for N = 233, Po = 0.29, A = 0.05 < A c , 
t = 300 < t c , and t = 1000 > t c . Modes are labeled as in (a), (d) Same, with relabeling as in (b). 
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Fig. 3 - Behavior of the deviations of the GAM from the initial value GAM(0) and of the largest 
discontinuity of the modulation function, (a) N = 144, A = 0.05, Po = 0.11. This velocity corresponds 
to the resonance 27rPo = ^2 9 /2. (b) Same N and A, but with Po = 0.1. The two values of Po are 
very close, but the behavior of the CM momentum is quite different since the second value is not 
a resonance. Nevertheless, the GAM and the largest discontinuity have the same kind of temporal 
evolution. This stresses the generality of the transition mechanism. 
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Fig. 4 - Behavior of t c i and t C 2 as a function of the number of particles N, for Po = 0.1, A = 0.05, 
and N = 144, 377, 610, 987, 1597. For all cases, t c2 > t cl . 



Indeed, writing the distortions x n as inverse Fourier transforms, with mode relabeling, gives: 

X n — S ' HmqC 9 (5) 
m 

From eqs. (Q) and (^) it follows immediately that: 

x n = f(nl) (6) 

This is true for the statics; taking into account the time dependence, we can generalize to the 
dynamics: a m — > a m (t), f(z) — > f(z, t). The equations of motion for the distorsions arel^ 



%n + Q — x-n+i + Xn—i — 2x n + A cos (nq + 2irx n + 2irQ) (7) 

From eqs. (||) and (Q), we derive the equation of motion for /: 

/ = /0 + l,t) + f(z -l,t)- 2/(z, t) - Q + A cos (2ttz + 2tt/(z, t) + 2ttQ(<)) (8) 

where we make the association x n ±x = z ± I. Let us now take the first of the two terms that 
make up the GAM in eq. (|^) and rewrite it in terms of the modulation function. If the mod- 
ulation function / is analytic, we can express the coefficients a nq as:a nq = dze~" l27TZ f(z, t) 
and using partial integration, we obtain: 

L = -i nqa- nq a nq = ^- f dz ^ f(z,t) (9) 

n 77 JO z 

We now compute the first time derivative of eq. (||) It can be shown that, due to periodic 
integration boundaries, only one term remains: 



L=^-\f dz ' cos (2tt(z + f(z, t) + Q(t)) (10) 
27T J dz 



Using / + z as a new integration variable and periodic integration boundaries, we get 



pi 

L = ~A / dzcos(27r(/(z,t) + z + Q(t)) (11) 



( 1 )In eq. (2) of ref. H and eq. (4) of ref. M the Q term has been omitted. 
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From eq. (0), we can now obtain the expression for the equation of motion of the CM Q: 

Q = — cos i.nq + 2lTX n + 2ttQ) (12) 

n 

Under the assumption that we made of the analiticity of /, we can go from the summation to 
the integral and, using eq. (||), we obtain: 

Q = X J d<j) cos {2ir{z + f(z,t) + Q(t)) (13) 

Comparison between eqs. ([ll]) and ( |l3| ) gives us (after integration): 

L + —Q = const (14) 

2-7T 

i.e., we have derived the conservation of the GAM from the definition of the modulation 
function. 

The crucial assumption in the calculation is the analiticity of the modulation function, 
which in the static case defines the regime under the Aubry transition. Notice that this result 
is compatible with our observation that t C 2 > t c \: analiticity of the modulation function 
implies GAM conservation, but the contrary is not necessary. Lastly, we notice from fig. || 
that this mechanism is robust: while resonances affect in a significant way the CM velocity, 
the GAM and the modulation function retain the same kind of behavior. 

Order of the transitions and relation with the static case. - Both transitions we observe 
in the dynamical case are very sharp; this suggests that, in the limit of an infinite chain, 
they could be first-order. This would be in contrast with the static case, where it has been 
proven (see e.g. ||) that the largest discontinuity in the modulation functions undergoes a 
second order transition. Presently, we do not have the analytical tools to check this rigorously. 
We can however get some numerical indications: we fit the behavior of the GAM and of the 
largest discontinuity to an exponential function, and study the behavior of the exponent as a 
function of the size N. In fig. ||, we see that, in both cases, the value of the fitting exponent 
grows when N increases. This indicates that, for N — > oo, the exponent could diverge and 
the transition would become first-order. More work is needed to confirm these initial results. 
It is also possible to relate in a direct way the dynamical and static behavior. This is shown 
in fig. [6] for i c2 : when the substrate potential approaches the critical static value A c , t c2 goes 
to zero. This allows us to relate the statics and the dynamics of the FK model. 

Conclusions. We presented in this paper numerical and analytical evidence that, in 
the kinetic FK model, a region is found where dynamical transitions take place. This region 
separates a floating IC from a pinned IC phase, and is therefore equivalent to the static 
Aubry transition (and, in fact, reducing to it in the limit of A — > A c ). Non- analiticity of the 
modulation function and breakdown of conservation of the GAM characterize these transitions. 
More work is needed to assess the order of these transitions and the relationship between t c \ 
and t c 2, investigating for example the possibility that, for N — > oo, t c \ and t C 2 converge to a 
single value t c . 
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Fig. 5 - (a) Exponential fits for the growth of the largest discontinuity in the modulation function 
for various N at A = 0.05, Po = 0.10. The initial values are determines by numerical accuracy, (b) 
Behavior of the fitting exponents shown in (a); the increase with N is apparent, (c) Exponential fits 
for the GAM. (d) Exponents of the GAM fit as function of N; here too, they increase with N. 
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Fig. 6 - Behavior of the transition time t C 2 as a function of the substrate potential A in the proximity 
of A c = 0.154.... Each line corresponds to a different initial velocity. t C 2 — > for A — > A c . 



